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Abstract 

We measured low energy cosmic-ray proton and helium spectra in the kinetic energy 
range 0.215 - 21.5 GeV/n at different solar activities during a period from 1997 to 
2002. The observations were carried out with the BESS spectrometer launched on 
a balloon at Lynn Lake, Canada. A calculation for the correction of secondary 
particle backgrounds from the overlying atmosphere was improved by using the 
measured spectra at small atmospheric depths ranging from 5 through 37 g/cm 2 . 
The uncertainties including statistical and systematic errors of the obtained spectra 
at the top of atmosphere are 5-7% for protons and 6-9% for helium nuclei in the 
energy range 0.5-5 GeV/n. 

Key words: cosmic-ray proton, cosmic-ray helium, atmospheric proton, solar 
modulation, superconducting spectrometer 
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1 Introduction 



Cosmic rays, charged particles from space, enter the atmosphere at a rate of 
several thousands per square meter per second. Among the cosmic-ray par- 
ticles, protons and helium nuclei are dominant components. Thus the energy 
spectra and absolute fluxes of these particles constitute fundamental data for 
studying cosmic-ray phenomena. The interstellar spectra carry important in- 
formation on the origins and propagation of cosmic rays in the Galaxy. At low 
energies, observable spectra in the heliosphere are deformed in a manner that 
depends on the solar activity. 

Cosmic-ray spectra observed at balloon altitudes suffer from a large back- 
ground of secondary particles produced inside the residual atmosphere. There- 
fore, in order to obtain the interstellar spectra and to understand the solar 
modulation, precise estimation of the atmospheric effects is important. The 
secondary components of protons are estimated by solving coupled transport 
equations (see for example Papini et al. [1]) or by performing Monte Carlo 
simulations assuming elementary nuclear interaction models between cosmic 
rays and atmospheric nuclei. The estimated secondary components can be 
checked and modified by comparing them with observed spectra at various at- 
mospheric depths. For particles below the geomagnetic cut-off rigidity, atmo- 
spheric cosmic-ray spectra consist of purely secondary particles, and therefore, 
cosmic-ray measurements at high geomagnetic cut-off rigidity at various alti- 
tudes are useful to improve the accuracy of the estimation of the atmospheric 
secondary particle contribution. 

We report here 0.215 to 21.5 GeV/n cosmic-ray proton and helium spectra 
measured by the BESS spectrometer for five balloon flights from Lynn Lake, 
Canada during the period 1997 through 2002. We also report on a realistic 
estimation of the atmospheric proton spectrum, which was tuned to be con- 
sistent with the proton spectra observed by the BESS spectrometer at small 
atmospheric depths. 

The low energy cosmic-ray spectra precisely measured at different solar ac- 
tivities are important to study the effect of solar modulation itself and also 
the cosmic phenomena observed at different solar activities. We already have 
much precise spectral data measured at balloon altitude or in space [2,3]. This 
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work, however, uniquely provides systematic measurements performed with 
basically the same spectrometer and at the same geographical condition. 



2 Spectrometer 

The detector for Balloon-borne Experiment with a Superconducting Spec- 
trometer (BESS) is a high- resolution rigidity (R = Pc/Ze) spectrometer with 
a large geometrical acceptance. It was designed [4,5] and developed [6,7,8,9] 
as an omni-purpose spectrometer to perform precise flux measurements of 
various components of cosmic rays [10,11,12,13,14], as well as highly sensi- 
tive searches for rare cosmic rays [15,16,17,18,19,20]. Fig. 1 shows a schematic 
cross-sectional view of the spectrometer in its 1997 configuration. The large 
acceptance is realized with a cylindrical structure and makes it possible to 
perform precise flux measurements with high statistics and small systematic 
errors. 

The spectrometer consists of a superconducting solenoid coil [5], a JET type 
drift chamber (JET), two inner drift chambers (IDCs), Time-of-Flight (TOF) 
plastic scintillation hodoscopes [8], and an aerogel Cerenkov counter [9]. A 
uniform magnetic field of 1 Tesla is produced by the thin superconducting 
coil which allows particles to pass through with small interaction probability. 
The magnetic field variation is less than 2.5% along a typical trajectory. The 
magnetic field region is filled with the central tracking detectors. Tracking 
of incident charged particles is performed by a circular fitting [21] of up to 
28 hit points, each with a spatial resolution less than 200 /zm, resulting in a 
rigidity resolution of 0.5% at 1 GV. The continuous, redundant, and three- 
dimensional tracking makes it possible to recognize background events with 
interaction or scattering. The upper and lower TOF counters measure the ve- 
locity (P) with 1/(5 resolution of 1.4%, and provide two independent ionization 
energy loss (dE/dx) measurements in the scintillators. Particle identification 
by mass is performed with these measurements(i?, f3 and dE/dx). Furthermore 
a threshold-type Cerenkov counter with a silica-aerogel radiator is installed 
below the upper TOF hodoscope. The radiator, with a refractive index of 1.03 
in 1997, and 1.02 in 1998 and later, was selected to veto e~ / ' jjT backgrounds 
for antiprotons up to 3.6 GeV in 1997, and 4.2 GeV in 1998 and later [9]. 
In the 2002 balloon-flight experiment, the observation was made with an up- 
graded detector developed for precise measurement of high energy particles 
up to 500 GeV [12,22]. 

The data acquisition sequence is initiated by a first-level TOF trigger which 
is generated by a coincidence of signals in the upper and lower TOF ho- 
doscopes with a threshold of more than 1/3 of the pulse charge for minimum 
ionizing particles. The trigger efficiency was evaluated to be 99.4 ± 0.2% 
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with a secondary proton beam at the KEK 12 GeV proton synchrotron [23]. 
The second-level trigger, which utilizes the hit patterns of the scintillator ho- 
doscopes and the IDC for the rough rigidity determination, selects negatively 
charged particles preferentially to improve the statistics for the antiproton flux 
measurement [18]. To build an unbiased trigger sample, one of every 60 (30 
in 2000 and 10 in 2002) events is recorded. In 1997, 1998 and 1999, another 
threshold level in the TOF trigger was set to be at 2.5 times higher pulse 
than that from minimum ionizing particles and one of every 25 events was 
recorded to improve statistics of helium events. To improve statistics during 
balloon-ascending periods in 1999 and 2000, we prepared a "low-energy proton 
trigger" which is generated by an on-board computer with a sampling rate of 
one out of two events to select particles with rigidity calculated to be lower 
than a fixed threshold. In the offline analysis, the "low-energy proton trigger" 
events are used to determine proton flux below 1.0 GeV during the ascent. 

In addition, an auxiliary trigger is generated by a signal from the Cerenkov 
counter to record particles above threshold energy with sampling rates of 
1/1 in 1998, 2000 and 2002, and 1/2 (1/1, 1/3, 1/4 during some periods) 
in 1999. This trigger improves the statistics for high energy particles. The 
efficiency of the Cerenkov trigger is evaluated as the ratio of the Cerenkov- 
triggered events among the unbiased trigger sample. It is higher than 90% 
for relativistic particles (/3 — > 1). For the flux determination, the Cerenkov- 
triggered events are utilized above 6.31 GeV (10.0 GeV in 2002) for protons 
and 3.98 GeV/n (5.41 GeV/n in 2002) for helium nuclei. Below these energies, 
the TOF-triggered events are used. 



3 Data Analysis 

3.1 Balloon flight observations 

We carried out balloon flights to measure cosmic-ray spectra with the BESS 
spectrometer at Lynn Lake, Manitoba, Canada (56°5'N, 101°3'W), where the 
geomagnetic vertical cut-off rigidity is low, approximately 0.4 GV. The bal- 
loons reached a level float altitude of 37 km. The residual atmosphere above 
the spectrometer was typically around 5 g/cm 2 . We also collected data during 
ascent in 1999, 2000 and 2002. A summary of the BESS balloon flights is listed 
in Table 1. 

We report here primary proton and helium spectra measured at Lynn Lake 
in 1997, 1998, 1999, 2000 and 2002. For a precise estimation of atmospheric 
secondary particle contributions, proton fluxes measured during the ascent 
period in 1999 and 2000, and during a slow descent period in 2001 were also 
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studied. A detailed analysis of secondary particles observed at Ft. Sumner, 
New Mexico, USA, where the vertical cut off rigidity is 4.2 GV, in 2001 is 
reported in Ref. [11]. 

During the period from 1997 to 2002, solar activity changed from the minimum 
to post-maximum as identified with sunspot number and neutron monitor 
data [24] as shown in Fig. 2. The magnetic field polarity also reversed during 
this period. 



3.2 Data reduction 

The procedure of data analysis and derivation of the flux was almost the 
same as that of the previous balloon flight data [10,11,12,13,14], except for 
the correction of atmospheric effects. 

We selected data sets which were taken with a stable condition of the spec- 
trometer and flight altitude. The live data-taking time of the selected data 
sets were 29531 sec, 25901 sec, 31890 sec, 41053 sec, and 28866 sec, in 1997, 
1998, 1999, 2000, and 2002, respectively. Among the selected data sets, we first 
selected "non-interacted" events, i.e., those associated with particles passing 
through the detector without interaction. The non-interacted event was de- 
fined as an event which has only one isolated track, one or two hit-counters in 
each layer of the TOF hodoscopes, and proper dE/dx inside the upper TOF 
counters. In order to estimate efficiency of the non-interacted event selection, 
a Monte Carlo simulation with Geant3 [25] was performed. The probability 
that each particle was identified as a non-interacted event was evaluated by 
applying the same selection criteria to the Monte Carlo events as that applied 
to the observed data. The resultant efficiencies of non-interacted event selec- 
tion at 0.3, 1.0 and 10 GeV/n were 82.7±0.9%, 83.2±0.9% and 76.5±1.1%, 
respectively, for protons, and 78.2±1.6%, 72.2±1.6% and 66.4±0.8%, respec- 
tively, for helium nuclei. The systematic error was estimated by comparing 
the hit number distribution of the TOF counters. Based on an accelerator 
beam test of the BESS spectrometer [23], the systematic uncertainty of the 
efficiency of non-interacted event selection for protons below 1 GeV should be 
within 2%. 

The selected non-interacted events were required to pass through the fiducial 
volume defined by the central region of the JET chamber, TOF counters and 
the aerogel radiator. In the analysis of ascent data, the zenith angle (9) of 
incident particles was limited within cos# > 0.95 to obtain nearly vertical 
flux. In order to check the track reconstruction efficiency inside the tracking 
system, the recorded events were scanned randomly, and the track reconstruc- 
tion efficiency was evaluated to be 99.5±0.2%. It was also confirmed that the 
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rare interacted events with more than one track are fully eliminated by the 
non-interacted event selection. 

3.3 Particle identification 

For the non-interacted events fully contained in the fiducial region, particle 
identification was performed as follows. The charge of the particle is identified 
by the ionization loss measurement. Both d-E/drc signals from the upper and 
lower TOF scintillators were required to be proton-like. Then, particles with 
proton mass were selected by a l//3-band cut. The selection of protons with 
the d-E/dx-band cut and l//3-band cut are shown on top and bottom in Fig. 3, 
respectively. Since at fixed rigidity the 1/(3 distribution is well described by 
Gaussian distribution and a half-width of the 1/(5 selection band was set at 
3.89 a, the selection efficiency was very close to unity (99.99% for a pure 
Gaussian distribution). Particle identification for helium nuclei was performed 
in the same way. However, 3 He nuclei were included in 1/ /3-band cut and were 
counted as helium-like events. In conformity with previous experiments, all 
the helium-like events were treated as 4 He nuclei in the analysis such as the 
reconstruction from rigidity to kinetic energy, and the efficiency estimation 
with the Monte Carlo simulation. 

3.4 Quality cut 

In order to maintain the quality of the rigidity measurement, a track quality 
cut was applied to the track fitting parameters for high energy particles above 
10 GeV/n. A TOF quality cut was also applied to confirm consistency between 
the hit position inside the TOF scintillator and the extrapolated track from 
the JET chamber. 

3.5 Contamination estimation 
3. 5. 1 Protons 

Protons were clearly identified below 1.7 GV by the mass selection as shown 
in Fig. 3. Above 1.7 GV, however, light particles such as positrons and muons 
contaminate the l//5-band for the protons, and above 4 GV deuterons (-D's) 
start to contaminate it. According to a Monte Carlo simulation of atmospheric 
cosmic rays [26], the flux ratio of (e + + fi + )/p is smaller than a few percent 
above 1.7 GV and decreasing with higher energy at balloon altitude. The 
observed D/p ratio was 2% at 3 GV by counting the number of events in each 
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l//3-band cut. The ratio at higher energy is expected to decrease [27] due to 
the decrease of escape path lengths of primary cosmic rays in the Galaxy [28]. 
The amount of those contaminations was as small as the statistic errors of the 
observed proton flux and no subtraction was made. 

3.5.2 Helium nuclei 

Helium nuclei were clearly identified by using both upper and lower TOF 
dE/dx selections. Since the l//3-band cut includes both 3 He and 4 He, the 
helium flux obtained includes both 3 He and 4 He. In conformity with previous 
experiments, all doubly charged particles were analyzed as 4 He. 

3.6 Corrections to the flux obtained at TOI 

The numbers of protons and helium nuclei passing through the BESS spec- 
trometer during the observation were obtained after correcting the detection 
efficiency. Then absolute flux at the top of the instrument (TOI) was ob- 
tained by taking account of energy loss inside the detector, live time, and the 
geometrical acceptance. The energy of each particle at TOI was calculated by 
summing up the ionization energy losses inside the instrument by tracing back 
the event trajectory. The total live time of data-taking was measured precisely 
by counting 1 MHz clock pulses with a scaler system gated by a "ready" status 
that controls the first level trigger. The geometrical acceptance was calculated 
by using the simulation technique [29]. In the high rigidity region where the 
bending of the track by the magnetic field is small, the geometrical acceptance 
for protons and helium nuclei is 0.160 m 2 sr without any limit for the zenith 
angle for the data at the float altitude, while it is 0.055 m 2 sr by limiting the 
zenith angle as cos# > 0.95 for the data during the ascent period. The simple 
cylindrical shape and the uniform magnetic field facilitate the determination 
of precise geometrical acceptance. The error arising from uncertainty of the 
detector alignment was estimated to be 1%. 



4 Corrections of atmospheric secondary particles 

The flux at TOI is measured under the residual atmosphere of about 5 g/cm 2 , 
and consists of a primary component surviving without any nuclear interac- 
tions with air, and of a secondary component produced inside the overlying 
atmosphere. In order to obtain the flux at the top of the atmosphere (TOA), 
the secondary component must be estimated and subtracted. The secondary to 
primary flux ratio for protons at 0.2 GeV is as high as 1.0 and 0.3 during solar 



7 



maximum and solar minimum, respectively, and the ratio decreases down to 
less than 3% above 10 GeV. Therefore the precise estimation of the secondary 
component is essentially important to obtain low energy proton flux at TOA. 
In this work we have tuned the estimation of atmospheric protons by using 
the observed data at balloon altitude as described below. For helium nuclei, 
the secondary to primary flux ratio is less than 5% and we used basically the 
same estimation as our previous works [10,12]. 

The proton spectrum inside the atmosphere can be estimated by solving simul- 
taneous transport equations following Papini et al. [1]. The primary spectrum 
at TOA is determined in an iterative procedure so that the estimated spectrum 
at TOI agrees with the observed one. Fig. 4 shows the comparison between 
the observed and estimated proton spectra at several atmospheric depths ob- 
served during an ascent period at Lynn Lake in 2000, and during a descent 
period at Ft. Sumner, New Mexico, in 2001 [11]. Fig. 5 shows the same com- 
parison but as a function of atmospheric depth. At Lynn Lake, the observed 
spectra contain both primary and secondary components, but at Ft. Sumner, 
the spectra below the geomagnetic cut-off energy of 3.4 GeV should be com- 
posed of purely secondary components. The dashed lines in Figs. 4 and 5 show 
the estimated spectra assuming the same production spectrum of atmospheric 
secondary protons as Papini et al. [1]. At all the observed atmospheric depths, 
the estimated spectra are much less than the observed data. 

In the energy region between 0.1 and a few GeV, the dominant source of the 
atmospheric secondary protons is recoil protons produced by interactions be- 
tween primary cosmic rays and air nuclei. We modified the energy spectrum 
of the recoil proton production in an iterative procedure to reproduce the ob- 
served proton spectra at several atmospheric depths between 5 and 37 g/cm 2 . 
The modified production spectrum was consistent with the original work by 
Papini et al. below 0.2 GeV as shown in Fig. 6. The solid lines in Figs. 4 and 
5 show the estimated spectra after the modification. The agreement between 
the estimation and the observation was significantly improved. We also con- 
firmed the agreement between the estimated and observed proton flux during 
an ascent period in 1999. The relative contribution of the secondary compo- 
nent is smaller in 1999 than in 2000, which was during the solar maximum 
period. The modified estimation can reproduce observed data irrespective of 
the residual atmospheric depth or sites (Ft. Sumner and Lynn Lake). 

In order to obtain the flux at the top of the atmosphere (TOA), secondary pro- 
ton production and interaction loss of primary particles inside the atmosphere 
were estimated by solving simultaneous transport equations with the modified 
production spectrum. The systematic error of the atmospheric correction for 
proton is estimated as 12.7%, which is composed of the uncertainty of resid- 
ual air depth (9.0%) and the cross section of primary cosmic rays with air 
nuclei (8.9%). The systematic error on the resultant proton flux is estimated 
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by multiplying the atmospheric correction factor as 3.7% and 7.6% in 1997 
and 2000, respectively, at 0.2 GeV. The systematic error becomes smaller for 
higher energy and is 1.2% and 1.6% in 1997 and 2000, respectively, at 1 GeV. 
The systematic error of atmospheric correction is as small as 1% for helium 
flux since the correction factor is less than 5%. 



5 Results and Discussions 

5.1 Proton and Helium spectra at the top of atmosphere (TO A) 

We obtained proton and helium spectra at the top of atmosphere during pe- 
riods of various solar activity from 1997 to 2002 as shown in Fig. 7. The 
numerical data are summarized in Tables 2 and 3. The overall uncertainties 
of the proton spectrum including statistical and systematic errors are 7.8% 
and 12.3% in 1997 and 2000, respectively, at 0.2 GeV, and 7.5% and 6.2% 
in 1997 and 2000, respectively, at 20 GeV. The overall uncertainties of the 
helium spectrum are 6.3% and 10.3% in 1997 and 2000, respectively, at 0.2 
GeV/n, and 18.2% and 7.7% in 1997 and 2000, respectively, at 20 GeV/n. 
The difference of the uncertainties at 0.2 GeV/n is due to the difference of 
the solar activity. The large uncertainties of helium spectrum at 20 GeV/n in 
1997 are due to the limited statistics. As described in Section 2, the Cerenkov- 
trigger was introduced from 1998 and the statistics were improved. As written 
in Section 3.5.2, all doubly charged particles were analyzed as 4 He to obtain 
helium flux. However the flux data contains a significant amount of 3 He. The 
isotope ratio ( 3 He/ 4 He) should be 10 - 25 % depending on energy [30]. 

The proton spectra at small atmospheric depths measured during the ascent 
period in 1999 and 2000 are summarized in Tables 4 and 5 and in Tables 6 
and 7, respectively. These data combined with the descent data in 2001 [11] 
were used to confirm our atmospheric corrections. 

A large solar modulation effect was found both in the proton and helium 
spectra in 2000 as a sudden decrease. On the other hand, proton and helium 
flux in 2002 slightly increased from those in 2000. We observed that the change 
in the flux variation between 2000 and 2002 was smaller than that between 
1999 and 2000 around a period of the solar magnetic field reversal. 

Since we have used an improved method for making the atmospheric correc- 
tions, this analysis is improved and should replace that done in our previous 
papers BESS-1998 [10], -1999 and -2000 [19], and -2002 [12]. The resultant 
spectra, however, are consistent with our previous ones above 10 GeV/n where 
the atmospheric correction is less than 3%. 
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5.2 Force Field Approximation 



In the Force Field approximation [31], cosmic-ray spectra at various solar 
activities are described by introducing the interstellar (IS) proton spectrum 
and one "modulation parameter" <fi. In this model, solar wind is assumed to 
be spherically symmetrical and is described with a simplified diffusion co- 
efficient and a transport equation for the propagation of cosmic-ray charged 
particles. The Force Field approximation predicts that the cosmic-ray particle 
of charge Z loses an energy of Z(f>, as if the particle decelerated in the static 
field with potential of <fi. Although the model is too simple to describe real 
dynamic activity of the solar wind, it has been generally used for the analysis 
of cosmic-ray spectra subject to solar modulation. The model is convenient to 
describe the basic feature of modulated cosmic-ray spectra and to indicate a 
degree of the solar activity. 

The proton and helium spectra obtained under this approximation are given 
in Fig. 8 and Fig. 9, respectively, in comparison with our measured data. The 
IS proton spectrum was assumed to be described with Af3 Pl R~ P2 ) where, as 
usual, (3 is the velocity of the particle divided by speed of light, R is the rigidity, 
and A, Pi and Pi are the fitting parameters. The parameter for BESS-1998 
was estimated to be ~600 MV by Myers et al. [32]. Other curves and values 
of given in Fig. 8 were obtained by fitting the measured spectra using the 
common IS proton spectrum. The Force Field approximation fits relatively 
well to the spectra measured in 1997, 1998 and 1999, which are in the positive 
phase of the Sun's magnetic field polarity. However, small discrepancies can 
be seen in 2000 and 2002, which are in the negative phase. According to recent 
works [33,34], the drift pattern of charged particles coming into the heliosphere 
varies with the sign of the polarity of heliospheric magnetic field. This feature 
cannot be treated in the Force Field approximation. Furthermore, the amount 
of energy loss depends on the observed particle energy, especially for negative 
polarity phase [34]. The small discrepancy seen in 2000 and 2002 may come 
from an inadequacy of the assumption of the model that the energy loss is 
independent of energy. We need further improvement in the models before 
we can consistently and accurately use them to estimate the shape of the 
interstellar spectrum. 



6 Conclusion 

Low energy cosmic-ray proton and helium spectra have been measured in 
the kinetic energy range 0.215 - 21.5 GeV/n with the BESS balloon flights 
in northern Canada during a period of solar minimum, 1997, through post- 
maximum, 2002. The proton spectra at TO A in the five flight experiments were 
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obtained after correction of the atmospheric proton fluxes measured in 1999 
and 2000. The correction was further improved by using the pure atmospheric 
proton flux measured in BESS-2001 at Ft. Sumner, where the geomagnetic 
cut-off energy is 3.4 GeV. The helium spectra at TOA were obtained in the 
same way. The overall uncertainties including statistical and systematic errors 
of the obtained spectra are 5-7% for protons and 6-9% for helium nuclei in the 
energy range 0.5 - 5 GeV/n. The maximum uncertainties are 12% for protons 
and 18% for helium nuclei in the whole energy range of the measurements, 0.2 
- 20 GeV/n. 

Assuming an interstellar (IS) proton spectrum and the modulation parameter 
of about 600 MV, the Force Field approximation can reproduce the measured 
spectra obtained from the BESS-1998 experiment, which are in the positive 
phase of the Sun's magnetic field polarity. Further corrections or model im- 
provements, however, would be required to reproduce the spectra during the 
negative polarity phase. The low energy cosmic-ray proton and helium spectra 
and the solar modulation effects have been much better understood based on 
the measurements with the BESS spectrometer in the solar minimum through 
post-maximum period. Those precisely measured data are important not only 
to study the effect of solar modulation but also to make a precise analysis of 
the cosmic phenomena observed at different solar activities. 
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z/tn J uly 


zytn July 


11-i-V\ A HIT 

ntn Aug. 


rutn Aug. 


z4tn aept. 


ith Aug. 


Float time [hr] 


on ^ 


oo n 


Q/1 ^ 


44. D 


1 /I Q* 

14. o 


1 P. K 
10. 


Observation time [hr] 


18.3 


20.0 


31.3 


32.5 


11.6 


13.9 


(during ascent) 






(2.8) 


(2.5) 


(2.4) 


(3.9) 


Atmos. Depth [g/cm 2 ] 


5.0 


5.5 


4.3 


4.3 


4.5-28 


4.6 


Reference 


[18] 


[10] 


[19] 


[19] 


[11] 


[12] 



* Including 11.8 hours of slow descent period 



Or 



Energy range 
(GeV nucleon -1 ) 



Flux± AFlux sta ± AFlux sys 
(m — 2 sr —1 s —1 GcV — 1 nucleon) 



BESS-1997 
[1997] 



; of BESS flights [observed year] 

BESS-1999 BESS-2000 
[1999] [2000] 



BESS-2002 
[2002] 



0.215- 


0.251 


1 


83 


±0 


03 


± 





13 


X 


10 3 


1 


22 


±0 


02 


± 





09 


X 


10 3 


9. 


88 


±0. 


20 


± 





76 


X 


10 2 


1 


69 


±0. 


04 


± 





18 


X 


10 2 


2 


22 


±0 


05 


± 


0. 


17 


X 


10 2 


0.251- 


0.293 


1 


.91 


±0 


02 


± 


0. 


13 


X 


10 3 


1 


.29 


±0 


02 


± 





.09 


X 


10 3 


1. 


04 


±0 


02 


± 





.07 


X 


10 3 


1 


92 


±0. 


04 


± 


0. 


18 


X 


10 2 


2 


.55 


±0 


.05 


± 


0. 


17 


X 


10 2 


0.293- 


0.342 


1 


92 


±0 


02 


± 


0. 


12 


X 


10 3 


1 


.35 


±0 


02 


± 





09 


X 


10 3 


1. 


12 


±0. 


02 


± 





.07 


X 


10 3 


2 


.06 


±0. 


04 


± 


0. 


17 


X 


10 2 


2 


.81 


±0 


.05 


± 


0. 


17 


X 


10 2 


0.342- 


0.398 


1 


92 


±0 


02 


± 


0. 


12 


X 


10 3 


1 


.35 


±0 


02 


± 





.09 


X 


10 3 


1. 


10 


±0. 


02 


± 





.07 


X 


10 3 


2 


13 


±0. 


04 


± 


0. 


17 


X 


10 2 


3 


.04 


±0 


.05 


± 


0. 


16 


X 


10 2 


0.398- 


0.464 


1 


84 


±0 


02 


± 


0. 


11 


X 


10 3 


1 


.33 


±0 


02 


± 





.08 


X 


10 3 


1. 


09 


±0. 


02 


± 





.07 


X 


10 3 


2 


.32 


±0. 


04 


± 


0. 


17 


X 


10 2 


3 


21 


±0 


.05 


± 


0. 


16 


X 


10 2 


0.464- 


0.541 


1 


.75 


±0 


02 


± 


0. 


10 


X 


10 3 


1 


.27 


±0 


02 


± 





.08 


X 


10 3 


1. 


05 


±0. 


01 


± 





.06 


X 


10 3 


2 


33 


±0. 


04 


± 


0. 


16 


X 


10 2 


3 


.28 


±0 


04 


± 


0. 


15 


X 


10 2 


0.541- 


0.631 


1 


.66 


±0 


02 


± 


0. 


09 


X 


10 3 


1 


.24 


±0 


01 


± 





.07 


X 


10 3 


9. 


85 


±0. 


13 


± 





.57 


X 


10 2 


2 


.43 


±0. 


04 


± 


0. 


16 


X 


10 2 


3 


.51 


±0 


04 


± 


0. 


14 


X 


10 2 


0.631- 


0.736 


1 


.49 


±0 


01 


± 


0. 


08 


X 


10 3 


1 


.16 


±0 


01 


± 





.07 


X 


10 3 


9. 


33 


±0. 


12 


± 





52 


X 


10 2 


2 


46 


±0. 


04 


± 


0. 


16 


X 


10 2 


3 


.52 


±0 


04 


± 


0. 


13 


X 


10 2 


0.736- 


0.858 


1 


.37 


±0 


01 




0. 


07 


X 


10 3 


1 


.07 


±0 


01 


± 





.06 


X 


10 3 


8. 


57 


±0. 


10 


± 





47 


X 


10 2 


2 


.48 


±0. 


03 


± 


0. 


15 


X 


10 2 


3 


46 


±0 


04 


± 


0. 


12 


X 


10 2 


0.858- 


1.000 


1 


20 


±0 


01 


± 


0. 


06 


X 


10 3 


9 


.79 


±0 


10 


± 





.52 


X 


10 2 


8 


32 


±0. 


10 


± 





44 


X 


10 2 


2 


49 


±0. 


03 


± 


0. 


14 


X 


10 2 


3 


.37 


±0 


.03 


± 


0. 


11 


X 


10 2 


1.00- 


1.17 


1 


.07 


±0 


01 


± 


0. 


06 


X 


10 3 


8 


.67 


±0 


09 


± 





46 


X 


10 2 


7. 


24 


±0. 


08 


± 





.38 


X 


10 2 


2 


.39 


±0. 


03 


± 


0. 


13 


X 


10 2 


3 


.28 


±0 


.03 


± 


0. 


11 


X 


10 2 


1.17- 


1.36 


8 


.93 


±0 


.08 


± 


0. 


46 


X 


10 2 


7 


.54 


±0 


.08 


± 





40 


X 


10 2 


6. 


60 


±0. 


08 


± 





.35 


X 


10 2 


2 


.34 


±0. 


03 


± 


0. 


13 


X 


10 2 


3 


.05 


±0 


.03 


± 


0. 


10 


X 


10 2 


1.36- 


1.58 


7 


.49 


±0 


07 


± 


0. 


38 


X 


10 2 


6 


.45 


±0 


.07 


± 





.35 


X 


10 2 


5. 


58 


±0. 


06 


± 





.30 


X 


10 2 


2 


.12 


±0. 


02 


± 


0. 


12 


X 


10 2 


2 


.87 


±0 


.03 


± 


0. 


09 


X 


10 2 


1.58- 


1.85 


6 


.09 


±0 


.06 


± 


0. 


31 


X 


10 2 


5 


.44 


±0 


06 


± 





29 


X 


10 2 


4. 


91 


= 


06 


± 





26 


X 


10 2 


1 


.98 


±0. 


02 


± 


0. 


11 


X 


10 2 


2 


.64 


±0 


02 


± 


0. 


08 


X 


10 2 


1.85- 


2.15 


5 


01 


±0 


.05 


± 


0. 


25 


X 


10 2 


4 


.47 


±0 


.05 


± 





24 


X 


10 2 


3. 


93 


±0. 


05 


± 





21 


X 


10 2 


1 


77 


±0. 


02 


± 


0. 


10 


X 


10 2 


2 


36 


±0 


02 


± 


0. 


07 


X 


10 2 


2.15- 


2.51 


4 


.03 


±0 


04 


± 


0. 


20 


X 


10 2 


3 


.62 


±0 


04 


± 





20 


X 


10 2 


3. 


11 


±0. 


04 


± 





17 


X 


10 2 


1 


.55 


±0. 


02 


± 


0. 


09 


X 


10 2 


2 


.05 


±0 


02 


± 


0. 


06 


X 


10 2 


2.51- 


2.93 


3 


.14 


±0 


.03 


± 


0. 


16 


X 


10 2 


2 


.90 


±0 


.03 


± 





16 


X 


10 2 


2. 


61 


±0. 


03 


± 





15 


X 


10 2 


1 


.34 


±0. 


01 


± 


0. 


08 


X 


10 2 


1 


.72 


±0 


01 


± 


0. 


05 


X 


10 2 


2.93- 


3.42 


2 


.39 


±0 


.03 


± 


0. 


12 


X 


10 2 


2 


.24 


±0 


.03 


± 





13 


X 


10 2 


1. 


98 


±0. 


03 


± 





11 


X 


10 2 


1 


14 


±0. 


01 


± 


0. 


07 


X 


10 2 


1 


.43 


±0 


01 


± 


0. 


04 


X 


10 2 


3.42- 


3.98 


1 


.83 


±0 


02 


± 


0. 


09 


X 


10 2 


1 


.71 


±0 


02 


± 





10 


X 


10 2 


1 


52 


±0. 


02 


± 





.09 


X 


10 2 


9 


.55 


±0. 


11 


± 


0. 


55 


X 


10 


1 


.17 


±0 


01 


± 


0. 


03 


X 


10 2 


3.98- 


4.64 


1 


39 


±0 


02 


± 


0. 


07 


X 


10 2 


1 


.27 


±0 


02 


± 





.07 


X 


10 2 


1. 


15 


±0. 


02 


± 





07 


X 


10 2 


7 


.75 


±0. 


09 


± 


0. 


45 


X 


10 


9 


.34 


±0 


.09 


± 


0. 


27 


X 


10 


4.64- 


5.41 


1 


.00 


±0 


01 


± 


0. 


05 


X 


10 2 


9 


.95 


±0 


15 


± 





.58 


X 


10 


9. 


28 


±0. 


15 


± 





54 


X 


10 


6 


.24 


±0. 


07 


± 


0. 


37 


X 


10 


7 


.30 


±0 


.07 


± 


0. 


21 


X 


10 


5.41- 


6.31 


7 


.54 


±0 


11 


± 


0. 


37 


X 


10 


7 


.08 


±0 


12 


± 





42 


X 


10 


6 


68 


±0. 


12 


± 





.39 


X 


10 


4 


.97 


±0. 


06 


± 


0. 


29 


X 


10 


5 


.59 


±0 


.06 


± 


0. 


16 


X 


10 


6.31- 


7.36 


5 


.33 


±0 


13 


± 


0. 


26 


X 


10 


5 


.01 


±0 


02 


± 





.30 


X 


10 


4. 


82 


±0. 


03 


± 





29 


X 


10 


3 


.74 


±0. 


02 


± 


0. 


22 


X 


10 


4 


.30 


±0 


.05 


± 


0. 


12 


X 


10 


7.36- 


8.58 


3 


.65 


±0 


10 


± 


0. 


18 


X 


10 


3 


.45 


±0 


01 


± 





21 


X 


10 


3. 


34 


±0. 


02 


± 





20 


X 


10 


2 


.70 


±0. 


01 


± 


0. 


16 


X 


10 


3 


11 


±0 


04 


± 


0. 


09 


X 


10 


8.58- 


10.0 


2 


.73 


±0 


.08 


± 


0. 


13 


X 


10 


2 


.41 


±0 


01 


± 





15 


X 


10 


2. 


43 


±0. 


02 


± 





15 


X 


10 


2 


.04 


±0. 


01 


± 


0. 


12 


X 


10 


2 


.22 


±0 


.03 


± 


0. 


06 


X 


10 


10.0- 


11.7 


1 


.86 


±0 


.06 


± 


0. 


09 


X 


10 


1 


.65 


±0 


01 


± 





10 


X 


10 


1. 


60 


±0. 


01 


± 





10 


X 


10 


1 


.42 


±0. 


01 


± 


0. 


09 


X 


10 


1 


.55 


±0 


01 


± 


0. 


05 


X 


10 


11.7- 


13.6 


1 


.25 


±0 


.05 


± 


0. 


06 


X 


10 


1 


.20 


±0 


01 


± 





.07 


X 


10 


1 


16 


±0. 


01 


± 





.07 


X 


10 


1 


.06 


±0. 


01 


± 


0. 


06 


X 


10 


1 


09 


±0 


01 


± 


0. 


04 


X 


10 


13.6- 


15.8 


8 


63 


±0 


.37 


± 


0. 


42 






7 


.96 


±0 


.05 


± 





49 






7. 


88 


±0. 


07 


± 





48 






7 


.22 


±0. 


04 


± 


0. 


44 






7 


.49 


±0 


04 


± 


0. 


27 






15.8- 


18.5 


5 


.71 


±0 


28 


± 


0. 


28 






5 


.61 


±0 


04 


± 





34 






5. 


52 


±0. 


06 


± 





34 






5 


15 


±0. 


03 


± 


0. 


31 






5 


.21 


±0 


.03 


± 


0. 


19 






18.5- 


21.5 


3 


.68 


±0 


21 


± 


0. 


18 






3 


.91 


±0 


03 


± 





24 






3. 


79 


±0 


05 


± 





23 






3 


.60 


±0 


03 


± 





22 






3 


.53 


±0 


03 


± 





14 







<<! to 

XS 

CD 
O 



o 

i-b 

Xi 

O 




o 

o 
xs 



CD 

B 
o 
t» 
xs 
ty 

CD 

1-1 

o 



CT> 



Energy range 
(GeV nucleon~ 1 ) 



Flux± AFlux 3 t a ± AFlux ays 
(m - 2 sr - ^s - 1 GeV - 1 nucleon) 



BESS-1997 
[1997] 



Name of BESS flights [observed 
BESS-1998 BESS-1999 
[1998] [1999] 



BESS-2000 
[2000] 



BESS-2002 
[2002] 



cr 



CD 

■ 'I >~ 



CO 



0.215- 


0.251 


2 


83 


±0 


08 


± 





16 X 


10 2 


2 


21 


±0 


07 


± 





12 x 


10 2 


1 


.78 


±0 


06 


± 





10 X 


10 2 


4 


78 


±0 


34 


± 


0. 


27 x 


10 


7 


66 


±0 


37 


± 


0. 


55 


X 


10 




0.251- 


0.293 


2 


.85 


±0 


.07 


± 


0. 


15 X 


10 2 


2 


.27 


±0. 


07 


± 





12 x 


10 2 


1 


.85 


±0 


06 


± 





10 X 


10 2 


5 


41 


±0 


.34 


± 


0. 


30 X 


10 


8 


.31 


±0. 


35 


± 


0. 


59 


x 


10 




0.293- 


0.342 


2 


90 


±0 


.07 


± 


0. 


15 X 


10 2 


2 


.29 


±0. 


06 


± 


0. 


12 x 


10 2 


1 


.90 


±0 


06 


± 


0. 


10 X 


10 2 


6 


32 


±0 


34 


± 


0. 


34 x 


10 


8 


54 


±0. 


33 


± 


0. 


61 


X 


10 




0.342- 


0.398 


2 


.75 


±0 


.06 


± 


0. 


15 X 


10 2 


2 


.16 


±0. 


06 


± 


0. 


12 x 


10 2 


1 


.89 


±0. 


.05 


± 


0. 


10 X 


10 2 


5 


.57 


±0. 


.30 


± 


0. 


30 x 


10 


7 


.55 


±0 


29 


± 


0. 


54 


X 


10 




0.398- 


0.464 


2 


.39 


±0 


.05 


± 


0. 


13 X 


10 2 


2 


.18 


±0. 


05 


± 


0. 


12 x 


10 2 


1 


.70 


±0. 


05 


± 


0. 


09 X 


10 2 


6 


02 


±0. 


.29 


± 


0. 


32 x 


10 


8 


.31 


±0. 


28 


± 





59 


X 


10 




0.464- 


0.541 


2 


.21 


±0. 


05 


± 


0. 


12 x 


10 2 


1 


.95 


±0. 


05 


± 


0. 


10 x 10 2 


1 


.55 


±0. 


04 


± 


0. 


08 X 


10 2 


6 


.30 


±0. 


.28 


± 


0. 


34 x 


10 


8 


.44 


±0. 


26 


± 


0. 




X 


10 




0.541- 


0.631 


2 


.05 


±0 


04 


± 


0. 


11 X 


10 2 


1 


.70 


±0. 


04 


± 


0. 


09 x 


10 2 


1 


.48 


±0. 


04 


± 


0. 


08 X 


10 2 


5 


.87 


±0. 


25 


± 


0. 


32 x 


10 


8 


.01 


±0. 


24 


± 


0. 


57 


X 


10 




0.631- 


0.736 


1 


.72 


±0. 


04 


± 





09 x 


10 2 


1 


.53 


±0. 


04 


± 


0. 


08 x 


10 2 


1 


.27 


±0. 


03 


± 


0. 


07 x 


10 2 


5 


.52 


±0. 


.23 


± 





30 x 


10 


7 


.23 


±0. 


21 


± 





51 


X 


10 




0.736- 


0.858 


1 


.38 


±0. 


.03 


± 


0. 


07 x 


10 2 


1 


.27 


±0 


03 


± 


0. 


07 x 


10 2 


1 


15 


±0. 


.03 


± 


0. 


06 X 


10 2 


5 


08 


±0. 


20 


± 





27 x 


10 


6 


.56 


±0 


19 


± 


0. 


47 


X 


10 




0.858- 


1.000 


1 


.23 


±0 


.03 


± 


0. 


07 x 


10 2 


1 


.09 


±0. 


03 


± 


0. 


06 x 10 2 


9 


.83 


±0. 


25 


± 


0. 


53 X 


10 


4 


.97 


±0 


19 


± 


0. 


27 x 


10 


6 


.16 


±0. 


17 


± 


0. 


44 


X 


10 




1.00- 


1.17 


9 


.36 


±0 


22 


± 


0. 


51 X 


10 


8 


60 


±0. 


22 


± 


0. 


47 x 


10 


8 


.16 


±0. 


21 


± 


0. 


44 x 


10 


4 


.34 


±0 


16 


± 


0. 


24 x 


10 


5 


.32 


±0. 


15 


± 


0. 


38 


X 


10 




1.17- 


1.36 


7 


.74 


±0. 


.18 


± 


0. 


43 X 


10 


7 


47 


±0. 


19 


± 


0. 


41 x 


10 


6 


.21 


±0 


17 


± 


0. 


34 x 


10 


3 


.80 


±0 


14 


± 


0. 


21 x 


10 


4 


.86 


±0. 


13 


± 


0. 


34 


X 


10 




1.36- 


1.58 


6 


.23 


±0 


15 


± 


0. 


36 X 


10 


5 


.57 


±0. 


15 


± 





31 X 


10 


5 


.32 


±0 


15 


± 





30 X 


10 


3 


.21 


±0 


12 


± 


0. 


18 X 


10 


4 


.09 


±0. 


11 


± 


0. 


29 


X 


10 




1.58- 


1.85 


4 


.66 


±0 


12 


± 


0. 


27 x 


10 


4 


.40 


±0. 


13 


± 





25 X 


10 


4 


.24 


±0 


12 


± 





.24 x 


10 


2 


.70 


±0 


10 


± 


0. 


15 X 


10 


3 


.25 


±0. 


09 


± 


0. 


23 


X 


10 




1.85- 


2.15 


3 


.73 


±0 


10 


± 


0. 


22 x 


10 


3 


.58 


±0. 


11 


± 





21 x 


10 


3 


.37 


±0 


10 


± 





19 X 


10 


2 


22 


±0 


09 


± 


0. 


13 X 


10 


2 


65 


±0. 


08 


± 


0. 


19 


X 


10 




2.15- 


2.51 


2 


.93 


±0 


.08 


± 


0. 


18 x 


10 


2 


72 


±0. 


08 


± 





16 X 


10 


2 


.56 


±0 


.08 


± 





15 x 


10 


1 


.85 


±0 


.07 


± 


0. 


11 X 


10 


2 


.19 


±0. 


07 


± 


0. 


16 


X 


10 




2.51- 


2.93 


2 


.17 


±0 


.07 


± 


0. 


14 x 


10 


1 


.90 


±0. 


07 


± 





11 X 


10 


1 


.91 


±0 


.07 


± 





11 X 


10 


1. 


.51 


±0 


.06 


± 


0. 


09 X 


10 


1 
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Energy range 
(GeV nucleon~ 1 ) 
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Fig. 1. Cross-sectional view of the BESS detector in its 1997 configuration. 
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Fig. 2. Counting rate of Climax neutron monitor (7-day average) and sunspot 
ber together with the BESS flight time [24]. 
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Fig. 3. Proton selection by d-E/dx-band cut for the top and bottom TOF hodoscope 
(a-1 and a-2) and l//3-band cut (b). 
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Fig. 4. Proton energy spectra at several atmospheric depths. Obtained spectra mea- 
sured during the ascent in 2000 (open markers), and the descent at Ft. Sumner in 
2001 [11] (closed markers) are compared with calculated spectra before (dashed 
lines) and after (solid lines) modification of the recoil proton production energy 
spectrum. 
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Fig. 5. Proton flux as a function of atmospheric depth. Obtained fluxes measured 
during the ascent in 2000 (open markers), and the descent at Ft. Sumner in 2001 
(closed markers) are compared with calculated fluxes before (dashed lines) and after 
(solid lines) modification of the recoil proton production energy spectrum. 
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Kinetic energy (GeV) 

Fig. 6. The original (dashed) and modified (solid) energy spectra of secondary proton 
production. For the original spectrum, two dominant components of evaporation 
and recoil production from air nuclei are also shown. The data points are found in 
Ref. [1] 
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Fig. 7. Proton and helium fluxes at the top of the atmosphere from 1997 through 
2002. Note that helium fluxes are plotted after being divided by a factor of 10. 
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Fig. 8. Proton spectra derived under the Force Field approximation. The interstellar 
(IS) proton spectrum was estimated by assuming <^>~600 MV for BESS-1998. Other 
curves and the value of 4> were obtained by fitting the BESS data. 
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Fig. 9. Helium spectra derived under the Force Field approximation. The interstellar 
(IS) helium spectrum was estimated by assuming </>~600 MV for BESS-1998. Other 
curves and the value of (f> were obtained by fitting the BESS data. 
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